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ABSTRACT 


In this thesis an attempt has been made to predict the 
utlimate lateral load capacity and the flexural behaviour of a 
pile when a certain depth of soil along the length of the pile 
has liquefied. The analysis has been carried out assumming 
that either there is complete or partial loss of soil support 
due to soil liquefaction over the estimated zone and computer 
programs have been developed for the same. The obtained 
results have been presented in the form of non-dimensional 
design charts for both short and long piles for various 
conditions of pile head fixity and ground conditions. 

The study brings out both qualitatively and quantitatively 
the significant influence of the liquefied soil zone on the 
lateral load capacity and bending behaviour of piles. The 
study also reveals that the presence of non-liquef ied soil 
cover increases the load capacity and decreases the pile top 
deflection in comparison to the corresponding values when the 
liquefaction occurs starting from the ground surface. This 
highlights the importance of soil densif ication in the upper 
zone of the soil. 
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Bending Moment Coefficient 
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Diameter of Pile, m 
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surface, m 
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Drag Force , kN 
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1. INTRODUCTION 


Liquefaction of saturated sands during earthquake is the 
cause of much damage to Civil Engineering Structures. Several 
instances of such damages have been observed in Niigata 
earthquake and Alaska earthquake that occured in 1964. The 
depth upto which a soil is likely to liquefy is governed by the 
ground water conditions, soil type, the earthquake magnitude 
and its epicentral distance from the concerned site. The 
liquefied depth may be estimated from different semi-empirical 
predictive models based on simple theories and experiments 
conducted in the laboratory and insitu to predict the 
liquefaction potential of a site. Some of these models are due 
to Seed et al. (1971), Lee et al (1972), Shibata et al. (1972), 
Ishihara and Yashuda (1972) , Yoshmi and Kuwabara (1973) , Yoshmi 
and Oka (1975) , etc. 

Once the probable depth of liquefaction is found out it 
can be used as an input data to estimate the effect of this 
liquefied zone on the existing foundations. Depending on the 
risks of soil liquefaction at a site different remedial 
measures like vibrof loatation , dynamic compaction etc for 
improving the ground condition may be adopted. Compaction of 
submerged sand to a relative density of 70 to 80 percent is 
adequate enough to prevent liquefaction. Alternatively the load 
bearing piles may be designed and adopted considering the 
effect of soil liquefaction on the pile capacity and the 
flexural behaviour under lateral loads. However, economy would 
be the main consideration in making the proper choice. 



Very often the engineers are called upon to make a 
comparative assessment of the various options available for 
making a final decision. Ground improvement techniques and 
processes have been studied over the years and design charts 
are available in literature. Considerable information on these 
techniques are now available in standard text books (van Impe, 
1989; Hausman, 1990) and reference books (Balasubramaniam et 
al. 1985) . Even though there are many theories available to 
predict the liquefaction potential at a site apart from damage 
survey the literature pertaining to the prediction of the 
effect of the liquefied depth on the foundation performance is 
scanty (Miura and O'Rouke, 1991) . 

For the cases when there is no soil liquefaction the 
ultimate lateral load capacity of piles may be estimated by 
using the methods proposed by Hansen (1961) and Broms (1964 (a) 
and 1964 (b) ) , Roy (1970), Meyerhof et al. (1981). Some of 
these methods are now available in standard text books on 
foundation engineering (Tomlinson, 1980; Poulos and Davis, 
1980) . Viggiani (1991) investigated the effect of sliding soil 
mass caused by land slides on the lateral load capacity 
considering various failure modes of the pile. The method has 
the potential to be modified and used in predicting the lateral 
load capacity of pile in liquefied soil deposit. 

The flexural analysis of a laterally loaded pile is 
generally carried out by using any one of the following 
approaches : 

(i) Matlock and Reese approach (1960) using the concept of 
modulus of subgrade reaction. 

(ii) Poulos and Davis approach (1980) using the Mindlin 
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solution. 


(iii) Finite element method alongwith appropriate constitutive 
modelling e.g. Desai and Abel (1972), Chen and Baladi (1985). 

Details of most of these methods have been presented by 
Poulus and Davis (1980), Zienkiewicz (1975), Desai and 
Christian (1977) and as such these are not reviewed here. 

It is evident that even though the detrimental effect of 
the loss of support due to soil liquefaction over a certain 
portion of an embedded pile is well recognised, the attention 
of geotechnical research community has not been focussed to 
quantify them. As such in this thesis an effort has been made 
to study the effect of liquefied soil zone on the ultimate 
lateral load capacity and flexural behaviour of a laterally 
loaded pile. 
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2. STATEMENT OF THE PROBLEM 


As shown in figure 2.1, a pile of length 1 and diameter d 
is embedded in a homogeneous saturated sand deposit. 
Considering the worst condition the water table has been 
assummed to be at the ground level. The depth of the soil 
cover which has not liquefied is Over the depth ± 2 , soil 

is in a liquefied state and the pile has lost the support of 
the soil either completely or partially. In the lower portion 
of the pile over a length I 3 the soil is again non-liquef ied. 
The sandy soil deposit is having an average submerged unit 
weight of r and the effective angle of shearing resistance 0 '. 
The horizontal load at the top of the pile is acting at a 
height e from the ground surface. The yield moment of the pile 
section is My. 

In this problem the possible movement of liquefied soil 
in the liquefied zone has also been considered. It has been 
assumed that the liquefied soil flows with a maximum velocity 
of Vq m/s at the top end of the liquefied zone as shown in the 
figure and lineraly varies to zero at its lower end. This 
flowing mass of soil will cause a drag force on the pile. 
However, the velocity of the soil movement is likely to be 
small. 

Thus as shown in the figure, for given inputs of 1^^, I 2 , 
I 3 , e and soil and pile properties, the objective is to 
determine the ultimate lateral load capacity of the pile, and 
the flexural behaviour of the pile when it has completely lost 
support over the specified length due to soil liquefaction. 
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3. ANALYSIS 


3*1 Ultimate lateral Load capacity 

In estimating the lateral load capacity the intensity of 
earth pressure acting at a point in a cohesionless medium due 
to the lateral displacement of the pile has been calculated by 
extending the Brom's theory (1964) and assuming that earth 
pressure intensities at any depth in the non-liquef ied zone can 
still be estimated by using the Brom's approach. 

As shown in the figure 2.1, the input data for the 
ultimate lateral load capacity of the pile are following: 

Ij^ = depth of upper non-liquef ied layer 

1 2 = depth of middle lic[uefied layer 

1 3 = depth of lower non-liquef ied layer 

Kp = coefficient of passive earth pressure of soil 
y = effective unit weight of the soil deposit 
P = density of the fluidised soil 
d = diameter of the pile 
My = yield moment of the pile. 

The various forces acting on a pile and their direction 
are shown in the figure 3.1. These are the forces when the 
failure is due to the yielding of the soil. 

a = soil resistance due to upper non-liquef ied layer 
= 1/2 * 3 y d Kp liXl^ 

= 1.5 y d Kp 1-^^ acting at 213/3 from the ground 

surface 

b = soil resistance (due to rectangular portion) for 
the lower non-liquef ied layer. 


6 


= 3 y d Kp + I2) . I3 at I3/3 from the upper end 


of 1- 


layer 


c = soil resistance due to triangular part 
= 1/2 * 3 y d Kp I3 . I3 

= 1.5yd Kp 1^3 at 2.I3/3 from the upper end. 

With reference to Fig. 2 . 1 , at any depth in the liquefied 
zone, the drag force FD can be computed by using the drag force 
equation for fluid mass as : 

FD = 1/2 p Cj3 . A 

where V = velocity of flow of the liquefied soil 
Cjj = drag coefficient 
A = Projected area of the pile 


Here 


^ = density of fluidised soil 

V = Vq X/ I2 


For an elementary length dx at distance x from the lower end 


d(FD) = 1/2 p Cjj.A 

= 1/2 P ( . xV 122 ) . Cp 


V. 


Cd • d 


d 
, 1 


dx 


FD 


1/2 


x2 . dx 


P^o_ •_Cd; d 

21^2 


P Vo d I2.CJJ 


If the location of the centroid of the force be at X2 
from the lower end of the liquefied zone then 


FD . X- 


r ^2 x2 

= d{FD).x = ( 1/2 P Vq — 

•’ Jo 


Cp.d.dx.x 
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P Vq^ .d.CQ.l2.X2 PVq^ ^ 2 ^ 

1/2 . 

6 122 4 

X2 = 0.75 I 2 

p V q d.l2.CQ 

Thus the drag force FD acts at 0.7 5 I 2 

6 

the bottom end of the liquefied zone. 

Now, using the above forces the estimation for the 
ultimate lateral capacity of pile can be made. The laterai 
load capacity for free head and fixed head piles have beer 
considered separately. 

3.1.1 Free Head Pile : To start with it is assumed that th< 
pile behaves. as a short pile and it fails due to the failure o; 
soil only. Here the maximum bending moment does not reach th€ 
yield moment of the pile section. 

The ultimate lateral load capacity of the pile for this 
case can be estimated by taking the moment of all the forces 
about the lower end of the pile (Fig. 3.1): 

ll I 3 I 3 

(e+1) + FD(0.75 I2 + I3) = a( + I2 + I3) + + 

3 2 3 

Hu = {a(l-0.67 13 ^)+ b. 13/2 +C.I3/3 -FD(0.75 I 2 +I 3 )} / (e+1) 

-(3.1.1) 

To check the validity of the assumption that the failure is not 
due to the development of plastic hinge in the pile, a check as 
follows is introduced in the analysis. 

If the maximum moment occurs in the 1^^ part and f is the 
depth of zero shear force below the ground surface then 
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or 


(3.1.2) 


= l.SydKpf^ 

f = 0.82 )V2 

ydKp 

Thus if f < 1]^ then the maximum moment occurs in the upper 
part, given by ( e + 2/3 f ) . 

Now, if < My then the pile is behaving as a short pile. 

On the other hand if » My then the pile behaves as a long 

pile. 

For the short pile case the ultimate lateral load capacity 
of pile is given by equation (3.1.1). 

For the long pile case the ultimate lateral resistance is 
given by the equations 

My = (e + 2/3 f) (3.1.3) 

H^=1.5rdKpf2 (3.1.4) 

If f calculated from equation (3.1.2) is such that f > 1], 
then the maximum moment will be in part I 3 . 

Let the maximum moment occur at a distance x from the 
lower end of the liquefied zone as shown in the figure 3 . 2 . 

The condition for maximum moment is that net shear force 
should be zero. Referring to Figure 3.2, 

+ FD = l.srdKpl]^^ + 3 y dKp (Ij^ + I 2 ) • X + l.sydKpX^ 

(here given from equation 3.1.1) 

The maximum bending moment occur ing at the zero shear depth is 
computed by taking moment of forces about that point. 

Mjjjax = (e+l 3 ^+l 2 +x) + FD (0.75 I 2 + x) - a (l]^/ 3 +l 2 +x) 

-3ydKp(l]^ + l 2 ).x.x /2 - 1.5 y dKpX^ .x/3 
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Mjjjax = Ha(e+li+l 2 +x) + FD (O. 75 I 2 +X) - a(lj^/3 + I 2 + x) 

-l.SrdKp ( 13 ^+ 12 ) x^ - 0.5 rdKpX^ 

Now again, if < My then the pile behaves as a short pile an 

its ultimate lateral capacity is given by equation (3.1.1). 

But if Mjjjax > My then the pile behaves as a long pile. 
For this case the ultimate lateral capacity will be given by 
the solution of two equations viz. 

My = (e + 13 ^+ 12 +x) + FD( 0 . 75 l 2 + x) - a(lj^/3 + I 2 + x) 

- 1.5 rdKp (l 3 ^+l 2 )x^ - 0.5ydKpx^ (3.1.5) 

and Hu + FD = 1 . 5 ydKpl 3^2 + 3 r dKp (l^^ + 12 )^ + 1-5 TdKpx^ 

- (3.1.6) 

The above two equations can be solved by an iterative method to 
give the ultimate lateral load capacity of the pile. This will 
also give the depth at which the yielding occurs. 

3.1.2 Fixed Head Pile : Assuming that the fixed head pile is 
behaving as a short pile, from fig. 3. 3 for the horizontal 
equilibrium of the pile 

Hu =a + b + c - FD 

where a,b,c have the same meaning as described earlier. 

Hu = 1.5 TdKplj^^ +3 r dKp( 13 ^+ 12 ) .I 3 + 1.57 dKpl 32 - FD .( 3 . 1 . 7 ) 

In fixed pile case, the maximum moment occurs at the fixed end. 
Thus taking moment of all the forces about the fixed end 
a(2/3 Ij^+e) + b(e+l-l 3 / 2 ) + c (e+l- 13 / 3 ) 

- FD(0.25l2 + li+e) (3.1.8) 

Now, if < My then the pile is behaving as a short pile and 

the ultimate lateral load capacity is given by equation 
(3.1.7) . 

) 0 



xr Mmax ^ "then the pile is behaving as an intermediate 
pile and the yield occurs first at the fixed end. 

Now from force equlibrium (Fig. 3. 3). 

F+Hu+FD=a+b+c 

where F is the reaction at the lower end of the pile and a,b,c 
and FD have the same meaning as described earlier. 

F = a + b + c - FD - Hu (3.1.9) 

Thus taking moment of all the forces about the fixed end of the 
pile. 

My = a(2/3 l^+e) + b(e+l-l 3 / 2 ) + 0 ( 6 + 1 - 13 / 3 ) 

- FD(0.25l2 + li + e) - F(l+e) 

From equation (3.1.8) the above equation can be written in a 
convenient form as 

^ = “max - F(l+e) (3.1.10) 

Now from equations (3.1.8), (3.1.9) and (3.1.10) the ultimate 

lateral capaity of the pile for the intermediate case can be 
computed . 

If the fixed head pile behaves as a long pile then another 
plastic hinge will form at some other point away from the fixed 
end. 

Now, if the second yield occurs in the I 3 part then taking 
moment about that point 

Hu (e + 2/3 f) = 2My (3.1.11) 

where f = 0.82 ( )-‘-/‘^ (from zero shear force ) (3.1.12) 

rdKp 

Thus the solution of these two equations (3.1.11) and 
(3.1.12) will give the value of ultimate lateral resistance for 
the long pile case. If the above calculated value of f > 
then the second yield may occur in the I 3 part. Again taking 



moment about that point which is at a distance x below the 
liquefied zone (Fig. 3.4) 

Hu (e+l]^+l 2 +x) - a(lj,/ 3 +l 2 +x) - 3 rdKp(l 3 ^+l 2 ) . xV2 

- 1 . 5 yd Kp X^/3 - FD(0.75 I 2 + X) = 2 My (3.1.13) 

The condition for maximum moment at that point gives 

Hu + FD=a + 3yd Kp(l]^+l 2 )x + 1.5 Y d Kp x^ (3.1.14) 

The solution of these two equations (3.1.13) and (3.1.14) gives 
the ultimate lateral capacity for the long pile case. 

Based on the expressions as deraved above a computer 
program has been developed to find the ultimate lateral load 
capacity of the pile when soil over a certain depth along the 
embedded lengh of the pile has liquefied. Using the program 
one can evaluate the lateral load capacity for both free and 
fixed head pile case and further find the failure mode of the 
pile i.e. whether it is failing as a short or as a long pile. 

For the purpose of providing non-dimensional design 
charts the drag force due to the possible movement of the 
liquefied soil has not been taken into account as the 
contribution due to drag force on the pile capacity is 
generally negligibly small. This can be quantitatively seen in 
the example as follows: 

For a free head short pile the expression for ultimate 
capacity is given by 

a(l-0.67lT) + b I3/2 + 1/3 C.I3 - FD(0.75 I 2 +I 3 ) 

. Hu = 

(e + 1 ) 

where a = 1.5 rd Kp l^^^ 

b = 3 y d Kp ( 13 ^+ 12 ) I 3 
c = 1.5 y d Kp I 32 



FD 


PV, 


dl2 Cjj 


Taking 


ll 

= 2 m 

^2 

= 2 m 

I3 

= 1 m 

1 

= 5 m 

r 

= 6 kN/m^ 

d 

= 0. 5 m 

e 

II 

0 

• 

01 

S 

= 3.0 

P 

= 1800 kg/m 


= 0.4 

a = 

1.5 X 6 X 0 

b = 

3 X 6 X 0.5 

c = 

1.5 X 6 X 0 


54 kN 


Fp = 1800 X Vq2 xO.5x2xO.4x 9.8 kN 


1000 


= 7.06 V, 


54(5-0.67x2)+ 108x1/2 + 2/3x13.5x1 - 7 . 06V^ (0 . 75x2+1) 


Hu = 


5.5 


260.6 - 17.6 
5.5 


= 47.4 


3.2 V, 


Taking = 0.5 m/s as soil movement is not likely to be very 
high. 


Hu = 47.4 - 3.2 X .5^ = 46.6 kN 

Thus we see that Hu value decreases from 47.4 to 46.6 only. 
Thus the drag force can be neglected without much error. 
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with a view to present non-dimensional design charts the 
equations as derived are non-dimensional ised as follows :- 
3.1.3 Non-Dimensional Form of the Equations 
Case I : Free Head Short Pile 

As previously derived the lateral load capacity for a free 
head short pile the expression is : 

H^(e + 1) = 1.5 rd Kp 1^2 (1 _ 2/3 Ij^) + 1.5 yd Kp (l^^ +12)13^ 

+ 0.5 yd Kp 13^ 

But I3 = 1 - 1^ - I2 
dividing both sides by ydKp 1^ 

1.5 (1]^/1)2 (1-0.67 li/l) + 1.5 (l^/l + I2/I) • 
_ (1 - li/l - l2/l)^ + 0.5 (1 - l^/l - 12/1)^ 

ydKpl2 (e/1 + 1) 


Case II - Fixed Head Short Pile 

H^ = 1.5 ydKp 1^^ + 3 ydKp ( 13 ^+ 12)13 + 1.5 yd KPI 32 
Since I 3 = 1 - Ij^ - I 2 . 

o 

= 1.5 (l3^/l)2 + 3 (li/l + I2/I) (1 - li/l - I2/I) 

ydKpl^ 

+ 1.5 (1 - l^/l - 12/1)^ 


Case III - Free Head Long Pile 

For the case of 1^^ = 0, from the previous equations 


H. 


u = 3 y d Kp 


lox + 1.5 y 


dKp 


x‘ 


and My = H^ (e + I 2 + x) - 1.5 ydKp l 2 X^ - 0.5 ydKpX^ 


H 


u 


ydKpl" 


3 (I 2 /I) 


(x/1) + 1.5 


(x/ 1 ) 



M H 

= ^ 2/1 + - 1*5 (I 2 /I) (x/l)2 

ydK 1^ rdKp 1 ^ ^ 

- 0.5 (x/l)^ 

Case IV - Fixed Head Long Pile 

For 12^=0 condition the equations are 
2 My = H^(e+l2+x) - 1.5ydKpl2x2 - O.SrdKpx^ 
and = 3 -ar dKp I 2 X + 1.5 r dKp x^ 

Thus 


2My . 

1 . = (e/1 +I2/I + x/1) 

ydKpi^ ydKpi^ 

j ^u 

and = 3 lo/l (x/1) + 

y dKpl^ 


- 1.5 l2/l(x/l)2 
-0.5 (x/1)^ 

1.5 (x/1) 2. 


Thus for the given parameters of I 2 /I ®/l the above 
equations can be solved to give the non dimensional lateral 
load capacity. 


3.2 Flexural Behaviour of the Pile 

Once the lateral load capacity of a pile is estimated for 
a given ground conditions the safe design load can be 
calculated by using a proper factor of safety and the flexural 
response of the pile can be analysed (Reese and 
Desai, 1977 ;Poulos and Davis , 1980) for its structural design. 
The adopted analysis procedure is briefly described as follows 

For deflection calcualtion the subgrade reaction analysis 
based on Winkler Soil model has been adopted ( Poulos and 
Davis, 1980) 

In the Winkler soil model the pressure p and deflection 
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« au Cl poim: are assumed to be related through a modulus oi 
subgrade reaction, which for a horizontal loading is denoted as 
kh- 

Thus p = kjjW . where k^ has the unit of force/length^ .The 

pile is assummed to act as a thin strip whose behaviour is 

d'^w 

governed by the beam equation El = -p.d. 

dx^ 

where E = modulus of elasticity of pile. 

I = moment of inertia of pile section. 

X = depth in soil, 
d = diameter or width of pile. 

Thus the governing equation of a laterally loaded pile is 
d^w 

El + kjjd . w = 0 . 

dx^ 

The solution to above equation can be most conventiently 
obtained by finite difference method ( Bowles, 1988 ) . The 
analytical solution are only available in convenient forms 
for the cases of constant and linear variation of kj^ along the 
pile. But in the case of cohesionless soil the value of kj^ is 
not constant and generally varies with depth. For homogeneous 
cohesionless soil deposits kj^ generally increases linearly with 
depth. To take care of the arbitrary variation of kj^ in the 
soil deposits it is convenient to correlate it with the SPT, N 
values. In this respect the variation of coefficient of 
subgrade reaction as proposed by the Japanese Highway Bridge 
Code (Yoshida and Hamada,1991) for such deposits has been 
adopted in the analysis as 

kj^ = 0.2 X 28 N (kgf/cm^) where d denotes the 

diameter of the pile in cm and N denotes the SPT value. 



set one tenth of that for non-liquefied layer. 

Now, the above basic differential equation can be written 
in finite difference form as given by Poulos and Davis, 1980 for 
a typical point i as (see figure 3.5). 


El Cw(i-2) - 4w(i-l) + 6 w(i) - 4w(i+l) + w(i+2) ] / 54 
+ = 0 

or, w(i-2) - 4w(i-l) + aj^ w(i) - 4w{i+l) + w(i+2) = 0 
where = 6 + ( kj^ 1^ d / El n^ ) 

where n = number of intervals along pile 

ki = modulus of subgrade reaction kj^ at ith node 
S = length of each subdivision along the pile 
1 = embedded length of the pile. 

The above equation can be applied from point 2 to n to 
give n -1 equations. 

Four further equations are obtained from the boundary 
conditions at the top and the tip of the pile. 

For free head pile the boundary condition at the top is 

d^w 

Shear = El = H 

dx3 

H 1^ 

or w(-2) + 2w(-l) - 2 W 2 + W 3 = (3.2.1) 

El n^ 


where H is the horizontal laod at top end of the pile, 
d^w 

Moment, El = H.e 

H. 

W(2) - 2w(l) + W(-1) = (3.2.2) 

El .n^ 

where e is the distance of the lateral load at the top of the 
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similarly for fixed head pile the rotation is zero at the 

top . 

Thus dw/dx = 0 

i.e. w (2) - w (-1) = 0 (3.2.3) 

Now, at the lower tip of the pile treating it to be a free end 

d^w 

Shear = El = 0 

dx^ 

i.e. - w(n~l) + 2w(n) - 2w(n+2) + w(n+3) = 0 (3.2.4) 

d^w 

and moment = El = 0 

dx^ 

i.e. w(n) - 2w(n+i) + w(n+2) = 0 (3.2.5) 

The final two equations are found from the overall equilibrium 
of the whole pile. 

For horizontal force equilibrium i.e. ^Fjj = 0 
As shown in the figure 3.5 the soil reaction at any node i 
is given by x A x w(i) 

= X ( S X d) X w(i) 

Thus Fjj ” 0 

H = ZQi 

or H = kj^ ( S X d) X w(i) ; i = 1, n+1 (3.2.6) 

or H = k-j^ (Sx d)/2 . w(l) + k 2 (8x d) w(2) + ^n+1 

(Sx d)/2 (W(n+i)) 

Now for moment equalibrium of the whole structure i.e. 

Xm = 0 

Taking moment of all the forces about the lower end of 


the pile. 



H (e+1) = Q^xL + Q 2 (L-S) + Q 3 (L-2S ) + Qn (S) 

+ Qn +1 • ( 0 ) 

it can also be written as 
gx d 

H(e+l)=ki ( ) l.w(l) + kjCSx d)(l-S).w(2) + k 3 (Sx d) (1-2S) 

2 

+ w(3) + +kjj (S X d) (S ) .Wj^ . (3.2.7) 

Thus equations (1), (2), (4), (5), ( 6 ), (7) together with n-1 

equations, apply for a free head pile and equations 1, 3, 4, 5, 
6 , 7 together with earlier (n-1) equations apply for fixed head 
pile. 

Thus, having n+5 equations in n+5 unknowns, these 
simultaneous equations can be solved by using Gauss elemination 
technique. The method is available in any standard textbook ( 
A. Ralston ,1965 ) on numerical methods of analysis. 

3.2.1 Non-Dimensional Form of Solution: 

The governing differential equation for pile can be 
written in a non-dimensional form and the flexural behaviour of 
the pile can be expressed in terms of non-dimensional 
parameters (Matlock and Reese, 1960) 

d^w 

El — + ki,.d.w = 0 

dx^ 

d'^w 

4 . _J w = 0 

dx^ El 

Putting K = kj^.d and defining a term relative stiffness as T = 
(EI/K) , which has a dimension of length. Thus the differential 
equation becomes (d^w/dx'* + T^ . w = 0 ) 

Defining non-dimensional terms as follows: 

Depth coefficient Z = x/T 



Maximum depth coefficient = L/T 

Soil Modulus function r = KT^EI 
Deflection Coefficient Ay = EI/H.T^ . w 
where H = Lateral Load at pile top. 

The differential equation can be written as 
d^Ay/dZ"^ + r . Ay = 0 


The corresponding boundary 
follows; 

At Top End of Pile 


d^w 

Shear : El 

dx-^ 

H.T^ 

Putting w = 

El 

d^Ay 

we get, — = l 

az^ 


d^w 

Moment : El 0 

dx2 

d^Ay 

we get, 

d z2 

At the bottom end of pile 


(3.2.8) 

can be written as 

= H 

Ay and Z = x/T 

(3.2.9) 

(3.2.10) 

= 0 


conditions 


Shear : d^Ay/dz^ = 0 (3.2.11) 

Moment : d^Ay/dz^ =0 (3.2.12) 

The equations (3.2.8), (3.2.9), (3.2.10), (3.2.11) and 

(3.2.12) can be written in the finite difference form and 
solved by Gauss elimination technique in the same manner as 
done previously. 

The non-dimensional bending moment coefficients Am is given 
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by 


Am 


/i2 

El d y , 

> mmf mm mm mmm * ^ mm mm mm mm 

H . T (1x2 

Thus, flexural behaviour of pile can be studied in terms 
of non-dimensional coefficients Ay and Am. 


2t 








— >j 4. — 




ZA 




cL 1^ 


FIG .3.4 
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4 . RESULTS AND DISCUSSION 


Based on the formulation prescribed in Section 3.1 for 
the ultimate lateral resistance of pile when certain depth of 
soil along the pile has liquefied, results have been provided 
in the forms of non-dimensional design charts for different 
cases of free head and fixed head piles. 

Representative results have also been presented to high- 
light the effect of soil liquefaction on the flexural behaviour 
of laterally loaded piles, based on the analysis as discussed 
in Section 3.2. Computations have been made using HP9000 
Computer System. The results and discussion pertaining to the 
lateral load capacity and flexural behaviour of piles are 
discussed separately as follows: 

4.1 Lateral Load Capacity : Figures 4.1 to 4.10 show the 
ultimate lateral load capacity of piles in terms of non- 
dimensional parameters which are defined as follows: 


H^/ y dKpl2 

M^/r dKp.l' 
I2/I 

ll/l 

e/1 


Non-dimensional lateral load capacity (NLLC) 

: Non-dimensional yield capacity of pile 
section (NYC) 

: Ratio of depth of liquefied zone to the 
total embedded length of the pile (NLD) 

: Ratio of depth of upper soil cover to the 
total embedded length of the pile (NCD) . 

: Overhang ratio with respect to the 
embedded length of the pile (OHR) 

In Figures 4.1, 4.2 and 4.3 the obtained results for the case 
of a free head short pile are presented. 

It can be seen from the figure 4.1 that when OHR and NLD 



are zero the well known Broiti ' s solution for lateral load 
capacity i.e. NLLC = 0.5 is obtained. The figures indicate a 
non-linear relationship between NLLC and NLD. Curves marked 1, 
2 and 3 in the figures 4.1, 4.2 and 4.3 show that as NLD 
increase NLLC decreases initially at an increasing rate, then 
at constant rate and finally at a decreasing rate. But curves 
marked 4 indicate continuous decrease in NLLC with NLD 
initially at a constant rate and then at a decreasing rate. 
These curves further indicate that in general for NLD value 
remaining constant as NCD increases NLLC decreases. This can 
be explained as follows: 

The increasing value of l^/l ratio (NCD) keeping I 2 /I 
ratio (NLD) constant indicates that the liquefaction is 
initiated increasingly at a greater depth and extends farther. 
As the greater part of soil resistance is mobilised at greater 
depths, this causes significant loss of soil resistance over 
the lower portion of the pile length resulting in lower values 
of pile capacity (NLLC) . But for the curve marked 
4 , corresponding to NCD equal to 0.4, the predicted values show 
the starting of a reverse trend beyond NLD equal to 0 . 4 . It can 
be seen that for any fixed 

value of NLD (0.6 > NLD > 0.4) as NCD increases upto 0.2 NLLC 
decreases but with further decrease in NCD, NLLC increases. 

The quantitative influence of l^/l (NCD) , I 2 /I (NLD) on 
lateral capacity (NLLC) is estimated in terms of a ratio R 
which is defined as follows: 

NLLC for given values of NCD and NLD 

R = 

NLLC with NLD = 0 

In Table 4.1 the values of this ratio (R) for different 
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NLD and NCD values are presented. To obtain the corresponding 
values of non-dimensional load capacity (NLLC) for OHR values 
of 0 , 0.1 and 0.2 the values presented in the table should be 
multiplied by 0.5, 0.45 and 0.42 respectively. The table 

brings out quantitatively the effect of the depth of soil 
cover (NCD) the depth of liquefied zone (NLD). on the lateral 
load capacity over the entire range of NLD from 0 to 1.0, and 
NCD from 0 to 0.4. 

In Figure 4.4 the obtained result for the case of fixed 
head short pile is presented. It can be seen from the figure 
that when I 2 /I (NLD) is zero, the Brom's solution for lateral 
load capacity for non-liquef ied case i.e. NLLC = 1.5 is 
obtained. Here also a non-linear relationship of NLLC and NLD 
is indicated. The NLLC values decrease at an increasing rate 
with increasing values of NLD. It is seen from the figure 
that for a given value of NLD, the NLLC decreases with the 
increasing value of NCD. The explanation of this behaviour is 
same as that of free head short pile case. 

Table 4.2 brings out the effect of soil liquefaction and 
soil cover on the lateral load capacity of piles 
quantitatively . 

In figures 4.5, 4.6 and 4.7 the variation of the non 
dimensional lateral load capacity (NLLC) with non -dimensional 
yield capacity (NYC) for the free head long pile case with 
different overhang ratio of 0, 0.1 and 0.2 are presented. These 
design charts have been drawn for the case when the cover ratio 
(NCD) is zero, and the liquefaction starts right from the 
ground surface. Thus the NLD=0, case in the present study 
(curve makred 1) is the same as that of Brom's (1964) non- 



liquefied case? The general trend of these curves indicate the 
following: 

- For a given OHR and NYC as NLD increases, NLLC decreases. 

- For a given OHR and NLD as NYC increases NLLC also increases 
at increasing rates. 

With the help of Table 4.3, for free head long piles a 
quantitative estimation of the influence of various parameters 
can be made. It is clear from the Table that for a increase in 
NLD from 0 to 0.6, there is a decrease in the load capacity of 
the pile ranges from 30 to 73 percent for different moment 
capacities of the pile. It is also seen from the table that 
for a variation in yield capacity from 0.013 to 0.42 (32 fold 
increase) the corresponding increase in load capacity is in the 
range of 10.3 to 19 times. 

In figures 4.8, 4.9 and 4.10 variation of NLLC with NYC for the 
fixed head long pile case, with different overhang ratios of 0, 
0.1 and 0.2 are shown. The general trend of the presented 
curves is similar to that of the free head long pile case viz . 

- For a given OHR and NYC, NLLC decreases with increasing NLD. 

- NLLC increases at an increasing rates with NYC for a given 
OHR and NLD. 

Table 4.4 gives a quantitative estimation of the 
influence of the various parameters on ultimate lateral load 
capacity for fixed head piles. It is seen from the table that 
for a 32 fold increase in NYC (from 0.013 to 0.42) the load 
capacity (NLLC) increase in the range of 8.13 to 29 times. For 
an increase in NLD from 0.0 to 0.6 the NYC decrease by about 
20-66 percent. 

An illustrative example for the studied range of values of NYC an 



OHR has been presented as follows to show the use of the develope 
non-dimensional design charts. 

Example : 

Given : 

Yield moment capacity of Pile section Mu =500 kNm 
Unit Weight of soil (effective) = 6 kN/m^ 

Cover depth, Ij^ = 0, 

Total embedded length of pile, 1 = 10 m. 

Passive earth pressure coefficient of soil, Kp=3 
Diameter of Pile, d = 0.5 m. 

Find : Ultimate Lateral Load Capacity of both free head and 
fixed head piles when 

Case A : (i) Depth of soil liquefaction I 2 = 4.0 m 

(ii) There is no liquefaction, I 2 = 0, and overhang 
length e = 0.0 and 1.0 m. 

Case B : Depth of soil liquefaction I 2 = 4.0 and the overhang 
length e = 0.7 m. 


Solution: 


From the given data : 


M 


u 


r dKpi- 


500 

= = 0.053 

10x0.5x3x10^ 


=0.4 

10 


Mu 


Free Head Pile : From Figure 4.5 for 


= 0.053 


ydKpl 


and 0.4 

1 
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when 


= 0.0 we get = 0.13 

ydKpl2 


====> Hu = 117 kN 


When e/1 = 0.1, from figure 4.6, 


0.11 


ydKpl" 


=> Hu = 99 kN 


For no liquefaction case i.e. when I 2 /I = 0*0 
For e/1 = 0.0, from figure 4.5 


= 0.25 


y dKpl' 


=> Hu = 225 kN 


For e/1 = 0.1 from Figure 4.6 


0.19 


ydKpl' 


=> H^ = 171 kN 


Fixed Head File From figure 4.8 (e/1 = 0.0) for I 2 /I = 


we get 


y dKpl^ 


0.25 


=> H^ = 225 kN 


From figure 4.9 (e/1 = 0.1) for I 2 /I = 


we get 


= 0.21 


ydKpl' 


=> Hy = 190 kN 


For no liquefaction case. 


For e/1 = 0.0 from figure 4.8 
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we get 


0.37 


we get 


Hu 

rdKpi^ 

=====> = 333 kN 

For e/1 = 0.1, from figure 4.9 

r dKpl^ 

=====> = 270 kN 


0.30 


Case B : 


For intermediate values of e/1 the value of Hu can be 
calculated by linear interpolation. e.g. for data given above to 
be the same if e/1 = 0.07 then Hu can be calculated as 
follows: 

For free head case 

Mu 

H^ = 117 kN, When e/1 = 0.0, I 2 /I = 0.4, = 0.053 

TdKpl^ 

Hy = 99 kN, when e/1 = 0.1, 
then for e/1 = 0.07, 

117-99 

H^ = 117 X 0.07 

0 . 1 - 0.0 

= 104.4 kN 

Similarly for fixed head pile 
H^ = 225 kN, e/1 =0.0 
H^ = 190 kN, e/1 =0.1 
For e/1 = 0.07, 


H^ = 225 


0 . 07 


200.5 kN 


225-190 

0 . 1 - 0.0 



4 . 2 Flexural Behaviour 


As already stated earlier the governing differential equation 
to find the flexural response of a pile in a non-dimensional 
form is: 

d^Ay 

+ r.Ay = 0 

dZ^ 

The boundary conditions for free head and fixed head pile case 
has already been given in section 3.2. The flexural behaviour 
of pile has been provided in terms of non-dimensional 
depth, deflection and bending moment coefficients for various 
cases of pile length, 

depth of liquefied zone and the soil cover over the liquefied 
soil. The non-dimensional terms used has been defined in 
section 3.2. The other terms used are as follows: 

Cover depth coefficient : Z-j^ = L^/T 
Liquefied depth coefficient : Z 2 = 

In the presented analysis over the depth Z-^, Z 2 and Z 3 
the original average SPT values when there was no liquefaction 
are taken as N^, k 2 N 2 ^ and k 3 N 2 respectively. When liquefaction 
occurs, SPT value is reduced and in the present analysis it has 
been taken as k 2 N 2^/10 over the liquefied depth as suggested 
(Yoshida and Hamada, 1991) . For computational purpose k 2 and k 3 
has been taken as 0.5 and 1.25 respectively. With these values 
the value of r in different zones would be 1.0, 0.05 and 1.25. 

The flexural behaviour of pile has been discussed in 
figures 4.11 to 4.26 for different values of cover depth Z-^, 
depth of liquefaction Z 2 , pile length Z^^. The different 
figures have been discussed separately as follows: 


4.11, 4.12 and 4.13 are for the free head pile 
case with Zj^ = 2.0. For this value of pile behaves as a 
short pile. These three figures together show the effect of 
soil cover Zj^ on the flexural behaviour of pile when the 
typical depth (Z^ + Z 2 ) upto which liquefaction is considered 
to be possible is constant and is equal to 0.8. Figure 4.11 
shows the deflection and bending moment variation along the 
length of the pile for the two cases when these is no 
liquefaction i.e. Z 2 = 0 (curve 2) and when there is 
liqpief action over Z 2 = 0.4 (curve 1) . Curve 3 and Curve 4 are 
similarly for bending moment variation. Figures 4.12 and 4.13 
show the variation of deflection and bending moment 
respectively along the length of the pile with varying from 
0.0 to 0.6. Based on the results presented in Figures 4.11, 
4.12 and 4.13, the ratio of the maximum deflection (Ay) and 
moment (Am) of a pile for different values of Z^ (soil cover) , 
Z 2 (non-dimensional liquefied depth) with the corresponding 
values for the non-liquef ied case are presented in Table 4.5. 
It is clear from the table that as the cover Zj^ decreases from 
0.8 down to 0.0 the relative values of the deflection 
coefficient Ay at the top of the pile increases by 5.4 times. 
It is seen from the table that with the decrease in the depth 
of non— liquefied soil cover the value of maximum bending moment 
coefficient Am, decreases upto a certain extent but when the 
cover is absent the Am value increases to its maximum value 
which is 2.6 times the value when there was no soil 
liquefaction at all. The above behaviour can be explained as 
follows: 
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For a given depth + Z 2 as the non-liquef ied 
cover depth decreases, the liquefied depth increases, the 
deflection of the pile also increases. This results in greater 
mobilization of soil resistance from the upper non— liquefied 
soil cover. This increased soil resistance in upper soil cover 
reduces the maximum bendinq moment in the pile at a lower 
section. But when the soil cover becomes very small, there is 
little resistance from the upper soil cover and there by 
having an insignificant influence on moment; so again the 
maximum bending moment coefficient value increases. 

Figures 4.14, 4.15 (a), 4.15 (b) and 4.16 are also for 
the free head pile case but with Z^j^ = 4.0 for which the pile 
may be treated as a long pile. These three figures show the 
effect of soil cover on the flexural behaviour of pile when 
the depth, (Z-j^+Z 2 ) is considered constant and is equal to 1.6. 
Figures 4.15 and 4.16 show the variation of deflection and 
bending moment along the length of the pile respectively. The 
results of figure 4.14, 4.15(a), 4.15(b) and 4.16 are presented 
in table 4.6. It is clear from the table that as the depth of 
cover decreases from 1.6 down to 0.0 the relative Ay value 
at the top of the pile increases by 5.95 times. The table shows 
that for this case also the maximum bending moment coefficient 
Am initially decreases with decreasing non-liquef ied soil 
cover. But when the soil cover becomes very small the Am value 
starts increasing. For no soil cover the maximum Am value is 
3.0 times than that for the case when there was no soil 
liquefaction. This behaviour can also be explained as above. 

Figures 4.17 and 4.18 are for the free head pile. Here 
maximum depth coefficient is equal to Z^ = 4.0 which indicates 



that the pile is a long pile. These two figures show the 
effect of depth of liquefaction Z 2 on the flexural behaviour of 
pile when depth of cover is constant and is equal to 0.4. 
Figures 4.17 and 4.18 show respectively the deflection and 
bending moment variation along the pile length. The results 
of figures 4.17 and 4.18 are shown in Table 4.7. It is clear 
from the table that as Z2 varies from 0.0 to 1.6 the relative 
value of Ay at the top of pile increases in the ratio of 1 to 
1.52. The table also shows that as the depth of liquefaction 
Z 2 increases for a constant cover depth the maximum bending 

moment coefficient Am decreases. It reduces to 0.62 times the 
initial value corresponding to no liquefaction. This 
behaviour can also be explained similarly as above. 

Figures 4.19, 4.20 and 4.21 are for the fixed head pile 
with Zjjj = 2.0. Thus it is a short pile. These figures show 
the effect of cover Zj^ on the flexural behaviour of pile when 
{Z-^-¥Z2) is constant and is equal to 0.8. The maximum deflection 
and bending moment occurs at the top of the pile. Table 4.8 
shows the relative variation of Ay and Am at the top of the 
pile. When the cover Z]^ varies form 0.8 down to 0.0 the 
relative Ay value at top of the pile increases by 1.95 times. 
The table shows that for given Z-^ + Z 2 as the non-liquef ied 
cover depth decreases the Am value first decreases and then 
starts increasing. The Am value for no soil cover case is 1.48 
times the Am value for a cover of Z^ = 0.8 i.e. when there is 
no liquefaction. This behaviour can be explained in a similar 
manner as that of the free head pile case. These variations 
are less as compared to the case of free head pile, thus 
showing the importance of the effect of fixity at the top of 



the pile for the fixed head case. 

Figures 4.22, 4.23 and 4.24 are for the fixed head pile 
with Zjjj = 4.0. These figures show the effect of on 
flexural behaviour of a long pile. Table 4.9 shows the 

relative variation in Ay and Am with As decreases from 

1.6 to 0.0 the Ay value increases by about 3.7 times. For this 
case also, the maximum Am value first decreases with 
decreasing non-1 iquef ied soil cover and finally starts 
increasing. Here the maximum Am value for no soil cover case 
is 2.2 times the Am value for the case of Z^ = 1*6 i.e. non- 
1 iquef ied case. 

Figures 4.25 and 4.26 show the effect of depth of 
liquefaction on the flexural behaviour of a fixed head long 
pile (i.e. Zj^ = 4.0). Here the constant cover provided is Z-^ = 
0.4. It is clear from the table 4.10 that as Z 2 increases 
from 0.0 to 1.6, the relative Ay values at the top increases by 
2 times and Am values increase by 1.07 times. The variation in 
maximum bending moment coefficient Am is very small as the 
depth of liquefaction increases from 0.0 to 1.6. 
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FREE HERD-SHORT PILE 

(NLD vs NLLC) 























FIXED HEnD-SHORT PILE 

(NLD vs NLLC) 
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FREE HERD -LONG PILE 

(NYC vs NLLC) 
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FREE HERD -LONG PILE 

(NYC vs NLLC) 
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FREE HERD -LONG PILE 

(NYC vs NLLC) 




FIXED HEnD-LONG PILE 

(NYC vs NLLC) 
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FIXED HEnD-LONG PILE 

(NYC vs NLLC) 
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FIXED HEnD-LONG PILE 

(NYC vs NLLC) 
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FREE HERD PILE 

{depth vs deflection, moment) 
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FREE HERD PILECZl Variable) 

(depth vs deflection) 
















FREE HERD PILECZl Variable) 

(depth vs moment) 















FREE HERD PILE 

{depth vs deflection, moment) 
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FREE HERD PILECZl Variable) 

(depth vs deflection) 
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FREE HERD PILECZl Variable) 

(depth vs deflection) 
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FREE HERD PILECZl Variable) 

(depth vs moment) 
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FIG 4-19 



FIXED HERD PILECZl Variable) 

(depth vs deflection) 
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FIG 4-20 






FIXED HERD PILECZl Variable) 

(depth vs moment) 
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FIXED HERD PILE 

(depth vs deflection, moment) 






FIXED HERD PILECZl Variable) 

(depth vs deflection) 
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FIG 4-23 
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(depth vs deflection) 
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FIG 4 2S 
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4.1 lioad Capaci1:y Variation for Free Head Short File 

( R = NLLC for given values of NCD and NLD / NLLC with NLD = 0) 


R 


NLD 

0.0 

0.2 

0.4 

0.6 

0.8 

0.9 

1. 

(NCD) 








o 

• 

o 

1 

0.90 

0.65 

0.36 

0.10 

0.04 

0 

H 

• 

O 

1 

0.82 

0.52 

0.25 

0.06 

0.04 


0.2 

1 

0.76 

0.46 

0.20 

0.12 

— 

- 

0.4 

1 

• 0.70 

0.46 

0.38 

.... 

—— 



Table 4.2 

: Load 

Capacity Variation for 

R 

Fixed Head 

Short Pile 

NLD 

0.0 

0.2 

0.4 

0.6 

0.8 

0.9 1.0 

(NCD) 







o 

• 

o 

1 

0.96 

0.83 

0.63 

0.36 

0.20 0 

0.1 

1 

0.92 

0.75 

0.53 

0.20 

0.03 

0.2 

1 

0.90 

0.68 

0.40 

0.053 

— - 

0.4 

1 

0.80 

0.53 

0.16 

— 

- 
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Tablo 4.3 : Variation of Non-’dimensional lateral load capacity 

for Free Head Long Pile 


NLLC 


NYC 

.013 

0HR= 

.42 

li 

1 



• 1 

013 

.42 

OHR=0 . 1 1 

% ' 

1 

.013 

.42 

1 

OHR=0 . 2 j 

(%) 1 

1 

(NLD) 

0.0 .075 

.85 

1 

(1033) 1 

.05 

.77 

1 

(1440) 1 

.04 

.68 

1 

(1600) 1 

0.2 

.05 

.80 

(1500) 1 

.03 

.70 

i 

(2233) 1 

1 

.02 

.65 

1 

(3150) f 

1 

0.4 

.03 

.70 

(2233) 1 

.02 

.60 

(2930) ' 

.016 

.55 

(3340) I 

f 

0.6 

.02 

.57 

(2750) 1 

1 

.015 

.53 

(3400) 1 

1 

1 

.015 

.48 

1 

(3100)J 

1 

1 

% 

-73 

-33 

1 

J 

1 

-70 

-32 

f 

1 

-62.5 

-30 

1 

f 


Table 4.4 : Variaton of Non-dimensional lateral load capacity 

for Fixed Head Long Pile 


NYC . 

013 

• 1 

O 1 
to 1 

1 

=0.0 1 
(%) 1 

.013 

.42 

0HR=0 . 1 1 

% 1 

1 

.013 

.42 

[ 

0HR=0 .2 1 
(%) i 
_Jl 

(NLD) 

0.0 . 

15 

1.37 

1 

(813)1 

1 

.10 

1.25 

1 

(1150) 1 

1 

.07 

1.12 

1 

(1500) 1 

0.2 . 

10 

1.30 

(1200) ] 

.07 

1.20 

(1614) ! 

1 

.05 

1.10 

(2100) 1 

1 

0.4 . 

06 

1.2 

1 

(1900)1 

.05 

1.10 

(2100) j 

.04 

1.00 

I 

(2400) ' 

0.6 . 

05 

1.05 

(2000)' 

1 

.04 

0.95 

1 

(2200) 1 

1 

t 

-I- 

.03 

0.90 

(2900) 1 

i 

1 

% 

-66 

-23 

1 

1 

1 

-60 

-24 

1 

1 

1 

-57 

-20 

1 

1 

1 
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Table 4.5 : Effect of Soil Cover on the Flexural Behaviour 

of Free Head Short Pile 
(Zj^+Z 2)=0.8, Zltt=2.0 


2l 

Ratio (Ay) top 

Ratio (Aia)jnax 

CO 

• 

o 

1.00 

1.00 

0.6 

1.06 

0.93 

0.2 

1.73 

0.88 

o 

• 

o 

1 

5.40 

2.60 


Table 4.6 : Effect of Soil Cover on the Flexural Behaviour 
of Free Head Long Pile 
(Zj+Zj) = 1.6, Zjj = 4.0 


Zl 

Ratio (Ay) top 

Ratio (Am)jj^ax 

1 

VO 

• 

H 

1.00 

1.00 

1.2 

1.03 

0.98 

0.8 

1.12 

0.85 

0.4 

1.45 

0.83 

0.0 

5.95 

3.00 
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TablA 4.7 X Effect of Liquefied Depth of Soil on the Flexural 

Behaviour of Free Head Long Pile 
“ 0.4, Zj^= 4.0 


^2 

Ratio (Ay) top 

Ratio 

0.0 

1.00 

1.00 

0.4 

1.18 

1.04 

0.8 

1.33 

0.97 

1.6 

1.52 

0.62 


Table 4.8 : Effect of Soil Cover on the Flexural Behaviour 
of Fixed Head Short Pile 

(Z^+Z2) “ 0.8/ Zjn “ 2.0 


Zl 

Ratio (Ay) top 

Ratio (Am) 

0.8 

1.00 

1.00 

0.6 

1.06 

0.96 

0.2 

1.44 

1.15 

0.0 

1.95 

1.48 



Tabltt 4.9 I Effect of Soil Cover on the Flexural Behaviour 


of Fixed Head Long Pile 
(Zj^+Zj) = 1.6, = 4.0 


Zi 

Ratio (Ay) top 

Ratio (Am)j,a 3 , 

1.6 

1.00 

1.00 

1.2 

1.12 

0.98 

o 

• 

CO 

1.25 

0.98 

0.4 

1.68 

1.05 

0.0 

3.70 

2.20 


Table 4.10 : Effect of Liquefied Soil Depth on the Flexur 
Behaviour of Fixed Head Long Pile 
ZjL = 0.4, Zjj = 4.0 

Z 2 Ratio (Ay) top Ratio 


o 

• 

o 

1.00 

1.00 

0.4 

1.22 

1.07 

CO 

• 

o 

1.50 

~1.07 

1.6 

2.00 

1 

• 

O 



5. CONCLUSION 


Based on the results and discussions presented in the 
Section 4 the following generalized conclusion can be drawn. 

(1) The lateral load capacity of piles decreases substantially 
with the increase in the liquefied soil depth. The decrease is 
more prominent when the liquefaction starts right from the 
ground surface. 

(2) Pile head deflections increase as the soil along the 
embedded length of the pile liquefies and is a function of 
various factors like the liquefied depth of the soil, length of 
the pile, soil cover, applied loads etc. 

(3) The beneficial effect of the ground improvement near the 
ground surface is significant. 

(4) The presented design charts and the computer program can he 
used fruitfully in determining the pile load capacity for short 
and long piles and finding the flexural response of a pil® 
under various conditions of soil liquefaction and the pile head 
fixity. 

(5) The presented design charts and the computer program are 
helpful in the quantitative analysis of the various factors 
affecting the pile capacity and the flexural behaviour under 
the various soil and pile conditions. 

(6) Providing fixed head pile is more beneficial under the 
conditions of soil liquefaction than free head pile. This 
observation is similar to the one made for the non-liquef i®*^ 
case. 
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6. SCOPE OP FUTURE STUDY 


(1) Analytical study for the determination of the lateral load 
capacity and flexural behaviour of a pile group under the 
various conditions of soil liquefaction and pile head fixity. 

(2) Experimental studies can be conducted with model piles to 
find the effect of soil liquefaction and soil movement on the 
pile behaviour. 
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APPENDIX 


fl-iSf f3nd^ the load capacity of a pile when the 

J I qu€ f ?rt 3 on orcured upto a certain depth of the pile. The 

Pile can b4 either © free head pile or a fixed heed pile. While 
er.ecutino tf'ie progra»n the relevant data is asked for. 

The variables !%sked for are pile di amet er , dept h where liquefac- 
tion tiartr i e It, depth of liquefaction 1£, total length of 
pile 1 . ec cen t r i c i t V e, unit weight of soil, passive earth pre- 
fsurf cofficient kp, yield noment of the pile my 
hroqra?rj starts from here 

reel d 1 1 , 1 & . 1 3 , 1 4 . 1 . e , g , k p , den , vo , cd , my , mma x 
read^, d }t,)?.l7,3,e,q,kp,den,vo,cd,my 

ft i Si'pt dt |vp » 1 t 3 f 

b *” 3 ♦ p t chM 3 t 1 ^ 1 3 h p 

c « I .S»g’*dtkp’M3't 13 
fd * dentvo>^vo^dtl£tcd/6 . 0 
f d f d/1 000 

wr i t e ( , f ) » j f free head pile then press 1 else 2' 
read’» , i t y pe 

program starts for free head piles 
assuming it to be a short pile 
if (it vpe . eq . 1 ) then 

hua ( ©T ( 1 <-0 . 67>H 1 )^b-t0 . 5^13ect:0 . 33^^ 1 3-f d^ ( 0 . 7S^ 1 2+1 3 ) )/(e+l ) 
to find maximum moment equate SF to zero 
assuminq maximum moment to be in the 11 part 
f * 0 . «82 + Eqrt ( hu/ ( d-t kp’^g ) 1 

i f ( f 1 e . 1 t ) then 

mmax « hu* ( e+0 . 67^f ) 
i f ( mma X . 1 e , my ) then 
hu *= hu 

write!*,*)' the pile is behaving as a short pile' 
wri t e (*,*)' t he maximum moment reached in kN is' 
wrlte(*,*) mmax 
else 

yt » my/( 1 . 5*g*d*kp ) 
f-0 

fs=f + 0.1 

f frrf 67*f ) 

if ( f I , It . y 1 ) then 
go to 50 
else 

write(*,*) 'pile behaving as a long pile' 
write(*,*)' and the pile has yielded' 
hu»^t 5*g*d*kp*f*f 
endi f 

endi f 

endi f 

maximum moment in 13 part 
if (f .gt .11) then 

c i = < a'-hu- fd)* 13*13 

X- ( sqrt ( ( b* 1 3 ) **2-4*c*c 1 ) -b* 1 3 ) /( 2*c ) 

mm a-'. =hu*te*l 1tl2*x) + fd*(0.75*12 + x >-a*(x + 12 + 0. 33*1 1 ) 

-0 . S* b* X * * 2/ 1 3 - ( 0 . 33* c* X **3 ) / ( 1 3* 1 3 ) 
i f ( mmax . 1 e my ) then 
hu«hu 

mmax = mmax *| 

write!*, ♦) ’pile behaving as a short pile' 
write!*,*) 'the maximum moment reached in kN-m is 
wri te < * , * ) mmax 
e 3 s f* 



( 

( 

c 

c 

<v. 

( 

( 

(>- 

L 

( 

L 

L 

C 

C. 

C 

o 


c yifldinp cace , b<?havl ng as a long plla 

y -•: (I 

60 y ■ y •» (I 0 t 

hu"< < 0 .33'* 1 1-*lc' + x ) + 0 . 5*b + x’*'*£/13 

* +0 33*r* y y»3/n 3t 2 J ) )/<e-t-l 1 •* 1£ + .'^ ) 

) = I a -hu~f d) * 1 3 * 1 3 

x?=(Egrt((btl3)*’*£-4*c*cl)-b*13)/(etc) 

5f(y gt.xl) than 

li= h> ! 

ur 2 1 1 > t pile ha« yielded The maximum' 

wr i t e ( moment has reached yield capacity' 
v.r i t e (♦,*)' t he depth of yield is x , x 
else 

0 0 to 6 0 
e n d i f 

d i f 
p T! d i f 

'the ultimate lateral capacity in kN is' 
wr i t e ( i ■♦= ) hu 
endlf 

lf(ltype.eq.2) then 

c here program starts for fixend end pile 

c 

c assuming shortplXe 

hu 1 f d 

inrnax=«a«<' ( e + 0 .6 7+11 )^b+ ( I'-O . 5+ 1 3 ) ( 1-0.33 + 13) 

♦ -fd+1e+n + 0 .2S+1S) 
lf(mmax.l©.my) then 
ha«hu1 

mma X X 

wr 1 1 eM , + ) ' t he fixed head pile behav'lng as short pile' 
wr i t e C + , + ) ' the maximum moment reached in kN-m is' 
wr 1 1 e t + , + ) mm ax 

wr 1 1 €( + .. + )' the lateral load capacity in kN Is' 

wr i t e( + , * 1 hu ' 

endlf 

1 f < mmax , gt . my ) then 
hU'S’hu 1 < my + mmax )/<©•♦ 1 ) 

c to check that maximum moment at some depth is in limit' 

• f « 0 . 8&+sqrt ( hu/( dakp=*g> ) 
i f ( f . 1 e . 1 I ) t e n 

mmax « hu+ 1 e+*0 . 6T + f ) 

1 f < mmax . 1 e . my I then 
hu » hu 

wrlte(»,+)' the pile is behaving as intermedate pill 
wri t e ( t , + ) ' the maximum moment reached in kN is' 
wrlteC+, ♦) mmax 
else 

yl «iEr+my/( 1 . 5 + g+d+kp ) 
f-0 

70 f^f^-O.I 

f ) 
if < f 1 . It .yl ) then 
go to TO 
else 

viritec*,*) 'pile behaving as a long pile' 
wr 3 te( + , + )'the pile has yielded at two points' 
hu=t .S+g’^d+kp+f + f 
endif 

endlf 

endi f 

c maximum moment in 13 part 

i f ( f et 1 1 ) then 

ct=<B-hu-fd)+13+l3 

( sqr t ( <t>‘M3)+ + E*^ + c + cl )-b+13)/(2+c) 



mill ax •^hU'»(e + 11 + l£-»'X)'*’fd»(0.75*l£+x)-a»<x + 12+0. 33* 11) 

-0 *b*x**£/13 -( 0 . 33*c*x**3)/( 13*13) 
i f ( mmax . 1 e . my ) than 

mmax -mmax 

wrlta(*,*) 'pile behaving as a short pile' 

Mrite(*,*) 'The maximum moment reached in kN-m is 

wriTeC*,*) mmax 

else 

yielding case j behaving as a long pile 

x = 0 

x==x+ 0 .01 

hu“(£*my+a* ( 0 . 33*1 1 +l£+x ) +0 . 5*b*x**2/13 
+ 0 . 33»c*x**3/< 13*13))/(ei-ll+l£+x) 
cl «( a-hu-fd)*13*13 
xl»(sqrt( (b*13)**£-4*c*cl )-b*13)/<2*c) 

1 f ( X . gt . X I ) then 
hu*hu 

wri t e <*,♦)' the pile has yielded at two points the 
wr i t e (*,*)' max i mum moment has reached yield value 
else 

go to 80 
endif 

endif 

endif 

write(*,*) 'the ultimate lateral capacity in kN Is' 
wriTe(*,*) hu 
endif 
endif 
stop 
end 

An example has been solved to show the use of the 
the program developed for the ultimate lateral load capacity 
for the liquefied case. The input data given is the following 
• The symbols used are the sa^he as that given in the program 
d » O.Sm, 11 « 3m, IS « 4m, 13 * 5m, 1 * 1£m, e » 0.5m, g « 
8kN/cu.m, Kp » 3, den » leOOkg/cu.m, Vo * 0.5m/s, Cd « 0.5 
my * SOOkNm. 

For this input for free head pile case the output comes as 
f ol lowi ng . 

'The pile has yielded. The maximum moment has reached 
yield capacity. The ultimate lateral capacity reached 
in kN is 174.41.' 
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c • 

c 

c 

c 

c 

c 

c 

c 


s 
1 0 

ciO 


c 


t 


30 


40 


SO 

c 


•> f d*?fl**ction ^nd bending inonient for lAterally load**d pile 
for d FREE HEAD head pile case. 

The variable's to be given are r1.»which is equal to 1/ tf~ 
which i-s the non-dimensional depth from the ground 'surface 
to where liquefaction starts, tS- which is the depth of liqu- 
efaction zm- is the total length of pile np- number of sub- 
divisions. 
parameter(nm3x=£00) 

double precision a ( £00 , £00 ) , p ( £00 , 1 ) , x ( 200 ) 
double precision d( £00 ) , ami £00 ) , av (£00 ) 
readt/rl ,zl , {£,/ffl,np 

do 10 i=l , np+1 
do 5 j=t , np+1 
a ( i , j )=0 
cont i nue 
coot inue 
do £0 i»l , np+1 
p ( t , 1 ) = 0 

continue 

np1 » i f i X ( z 1 ♦ np/zm+1 ) 
np£« i f i X ( ( z 1 + z£ ) *np/zm+1 > 
del*zm/np 

program 

r£»rl 

r3-l .£S*r1 

p( 1 , 1 )*del + + l 
a(l,1 )*(e+rl»del* + 4) 
a( 1 ,£)*-4 
a( 1 ,!)*£ 

*(£,1 )«-£ 

a{2,2)*(S+rl »del**4) 

a(2,3)*-4 

•»l£,4)»1 

do 30 i=3,npl 

a( i, i-£)=l 

a I i , i-1 )*-4 

a<i,i)=(6+r1»dele+4) 

a( i , i+1 )*-4 

a( i , i + c')»l 

continue 

do 40 i»npl +1 , npE 

a(i , i-£)=l 

a ( i , i -I )*-4 

aii, i )*(6+r£+del*«A) 

ai 1 , 1 + 1 >*=-4 

a( i , i+£)=1 

contltiue 

do 30 i«np£+l,np-l 
a i i , i -£ ) * I 
aii , 1-1 )=-4 
ai i , i )«( 6 + r3 + del**4 ) 
aii , i + 1 )»-4 
aii , i + £)«1 
cont inue 
*»*+»e+e*» 
a(np,np-L*)»l 
ainp, np-1 )»-4 
4 ( np , np ) « ( S+r3*de 1 **4 ) 

4i np , np+ 1 ) »-£ 

^ d. it; ± J* A 



« * in ( I ) * 0 

do 51 t?, n[> 

dm(i ) = (x(tf 1 )-2*x ( i ) + x( i-l ) 

SI con t I nu*» 

i»m ( n|7+ 1 ) »i) 

•4V ( 1 ) * 1 

do 5£ i *1 , np- 1 

,av( i ) = ( - X { i )+e*x ( i-1 )-£»x ( i + n + xCi+Sll/Cdel**!) 
Sc! con t i nu» 

av ( np ) =0 
av ( np+ 1 ) *0 
do 53 1 = 1, np'»- 1 
d(i)=(i-1 l+d-l 
53 cont i nu** 

do 60 i= 1 , np+1 
wrii:e(*,61) d(i>,am(i) 

6 0 ront 1 nut* 

61 f ormat ( £x , £f 9 . 3 ) 

stop 
end 

subroutine qauss ( a , p , x , n ) 

double precision a(£00,£00 ) . p(Si)0, 1 ) , x (SOO ) 
do £00 i = 1 , n, 
d o £ 0 0 j = i , n 

if ( i . ne . j ) then 

p< j, I )=p(j,1)-a(j,i)»p(i, })/a(i,i) 
do 100 k=n , 1 , - 1 

10 0 a ( j , k ) *a ( j , k ) - j ( i , k ) * a ( j , 1 ) /a ( i , i ) 

end i C 

£00 continue 

do 3 00 i = 1 , n 

300 X ( 1 )= p( i , I )/a( 1 , i ) 
reTurn 
end 



c 

* 

c 

c 

c 

c 

c 

c 


1 0 


c'O 


c 


c 


30 


40 


SO 

c 


•>f d**fJ»»ct ion jnd binding momi'r)i: for t iojoho pi., 

for A FIXEO HEAD h^-td pilt c«s* . 

Tim v«ir i ub I i*s to b* giv^n •r* r1, which is »qu 4 l to 1, z1“ 
whicli i-i the non-dimensional depth from the ground -surfHice 
to where I i quef act i on starts, which is the depth of liqu- 
efaction zin- is the total length of pile np- number of sub- 
divisions. 
paramet«r(ntndxs£00) 

double precision a ( £00 , £00 ) , p ( CO 0 , 1 ) , x ( £00 ) 
double precision d ( 800 ) , am ( £00 ) , av ( £00 ) 
read*,rl , z1, <£,zm,np 

do 10 i*1 , np+1 
do 5 j = 1 , np+1 
a ( i , j ) * 0 
cont inue 
cout inue 
do £0 i=l , np+1 
p ( t , 1 ) = 0 

con t i nue 

npl *i f i X ( i 1 +np/zm+l ) 
np£=i f ix ( ( z1 +z£)*np/zm+l ) 
del»zm/np 

program 

r£*r1/£«) 
r3»1 . £S*r 1 

p( 1 , 1 ) “d** I »♦! 
a ( 1 , 1 ) =( 6+r 1 *del *+4 ) 
a(l,2)»-0 
a(1,3)*£ 
eeeyee 
a(2,1)»-4 

4(2,8)*(7+r1 +del**4) 

a(2, 3)«-4 

a(£,4)=1 

do 30 i=3,np1 

a(i,i-2i=l 

a( i , i-1 )=-4 

a( i , 1 )*( 6 + rl *del*»4) 

a( i , t + 1 )=-4 

a ( i , i + 8 ) - I 

cont I nue 

do 40 i * np 1+ I, np£ 

a{i,i-£)=1 

a(i,i-l)=-4 

a { i , i ) » ( 6+r£*del * *4 ) 

a( i , i+ I )*-4 

a(i, i+£)=1 

continue 

do 50 i *-np£+ 1 , np-1 

a(i,i-8)=l 

a(i,i-1 )=“4 

a(i,i)*(e+r3»del**4) 

a( i , i+1 )=-4 . , 

a( i , i+£ )* 1 

continu*? 


.» ( np, np-£ ) » 1 
a(np,np-l )*-4 
a(np, np )«(S+r3+del**4) 
a(np,np+1 )»-£ 



^ « in ( t ) « 0 

do 51 t« £,np 

din(i) = (x(t + l )~S*x(i1 + x(i-l ))/(d*»l*»£) 

SI continue 

im ( np + 1 ) *0 
*v(1)«l 
dv ( £ ) »1 

do 5 £ i ■ “3 , 11 p ~ 1 

4v(t)»(-*K{ l-£)+e*x(i-1 )-2*x(i + 1 ) + x(l+2) )/(d*l**3) 
S£ con t 1 nu<» 

*v( np)»0 
iv ( np+1 ) =0 
do 53 i»1 , np-i-l 
d( i )«{ i-1 l^d.*! 

53 continue 

do 60 i“1 , np+1 
write(*,61) d(i)>am(i) 

60 continue 

61 format (£x,2f9 ,3) 
stop 

end 


100 

200 

300 


subroutine g&usst a, p , x j n ) 

double precision a < 200 , 200 ) , p ( 200 , t } , x (200 ) 
do 200 i-1 ,n ^ 
do 200 j-1 , n 

if ( i . ne . j ) then 

p(j, 1 )®p(j,1 )-a(j,i)*p(i,1 )/a(i,i) 
do 1 00 k=n, 1 , “1 

a(j,k)sa(j,k)-a(i,k)»a(j,i)/a(i,i) 
endif 
cent i nue 
do 300 i*l,n 
x( i )- p( 1 j 1 )/a{ i, i ) 
return 
end 



